Dislocations Jam At Any Density 
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Crystalline materials deform in an intermittent way via dislocation-slip avalanches. Below a criti- 
cal stress, the dislocations are jammed within their glide plane due to long-range elastic interactions 
and the material exhibits plastic response, while above this critical stress the dislocations are mobile 
(the unjammed phase) and the material flows. We use dislocation dynamics and scaling arguments 
in two dimensions to show that the critical stress grows with the square root of the dislocation 
density. Consequently, dislocations jam at any density, in contrast to granular materials, which only 
jam below a critical density. 
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When a crystalline material is sufficiently deformed, 
it undergoes irreversible, or plastic deformation. Tra- 
ditionally, plastic deformation of crystalline solids has 
been considered to be a smooth process in time, and ho- 
mogeneous in space, since fluctuations are expected to 
average out at sufficiently large spatial scales. On small 
spatial scales, however, intermittent motion of disloca- 
tions is observed, resulting in a pattern of deformation 
that is spatially inhomogeneous and intermittent in time. 
This behavior is associated with a coherent motion of the 
dislocations that releases stress by slip avalanches: se- 
quences of events with long-range correlations in space 
and in time. The slip avalanches span several orders of 
magnitude in size and the energy released is distributed 
according to a power law [THS]. 

Theoretical models, including discrete dislocation dy- 
namics models [U ITHTT]. continuum models [H fTTHl3] . 
phase field models [HI [M] and phase field crystal mod- 
els [TS] are able to reproduce many of the experimental 
findings and reveal scale invariant, power-law distributed 
phenomena that are indicative of a non-equilibrium criti- 
cal point [E]. The dislocation system is jammed below a 
critical value of the external stress. Applying a constant 
external stress above the critical (yield) stress allows the 
system to flow, and the dislocations are unjammed. It is 
important to stress that in the glide plane of the dislo- 
cations, there is effectively no external potential, so that 
the jamming is an emergent phenomenon. However, re- 
cent work has shown that the behavior of the transition 
appears to be in the universality class of the interface 
pinning-depinning transition 115j , as if there was an 
effective external potential induced by the collective in- 
teractions between the dislocations. 

In this Letter we study connections between the plas- 
tic yield point of systems with long range interactions, 
such as dislocation systems and the jamming transition 
of systems with short range interactions, such as sheared 
granular materials and molecular liquids [TTH^n] . When 
a system jams it undergoes a transition from a flowing 
state (analogous to a depinned phase) to a rigid state 
(analogous to a pinned phase). In contrast to the or- 
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FIG. 1: (color online) Proposed phase diagram for dislocation 
systems. Notice the absence of a jamming point. 

dered solid phase obtainable via crystallization, the solid 
phase reached via jamming remains disordered. Liu and 
Nagel Tl\, O'Hern et al. [18l[l9] and others [20] studied 
jamming of granular materials with short range interac- 
tions in simulations and experiments. They found that 
below a critical density these materials do not jam at any 
stress. This critical density is called the jamming point 
J of granular materials. In contrast, we show here that 
dislocations jam at any nonzero density, i.e. dislocations 
have no jamming point. The physical reason is that dis- 
locations have long range interactions that can lead to 
pinning for arbitrarily large distances between the dislo- 
cations. Fig. 1 sketches the putative jamming phase 
diagram (in the absence of screening) for dislocation- 
mediated plasticity. It is closely related to the jamming 
phase diagram of [5D] for granular materials, except for 
the absence of a jamming point J for dislocations. 
In the following we employ analytical calculations and 
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discrete dislocation dynamics simulations to study how 
the critical yield stress depends on the dislocation density 
p. Our analytical calculations verify and generalize the 
numerical findings. 

The Model:- We place N straight edge dislocations par- 
allel to the 2;-axis in a square box of side L. They are 
allowed to glide only along the shear direction (x-axis), 
while they can interact in the x and y directions. This 
simulates a single-slip system. Materials, like ice, with 
strong plastic anisotropy, deform by glide on a single 
plane [T]. In these systems dislocation climb is negli- 
gible due to high plastic anisotropy. The z direction has 
been shown to be irrelevant to scaling [l][7HT0] effectivelly 
rendering the problem two-dimensional. An edge dislo- 
cation with Burgers vector b — (6, 0) produces in the host 
medium an elastic shear stress at a distance r — {x,y), 
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where is the shear modulus and v is the Poisson ratio 
of the host medium [5T]. This is anisotropic in the (x,y)- 
plane and decays as Tint ~ ^/r at large r = \/ x"^ ^ y"^ . If 
an external shear stress Text is applied the overdamped 
equation of motion of a dislocation along the shear direc- 
tion is described by, 
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for z,j = l,...,iV where Xi is the x coordinate of the 
«th dislocation at point with Burgers vector hi , Vj with 
2 ^ i are the coordinates of the other iV— 1 dislocations, t 
is time and r\ is the effective viscosity in the host medium 
[UlTllH]. Here we have set the temperature to T = 0. For 
the computer simulations we have set the distance scale 
6=1 and the time scale = '?/(m/(27I'(1 — v)) — 1. 
To simulate bulk materials we employ periodic boundary 
conditions in both x and y directions. 

To treat the long-range character of the dislocation in- 
teraction, we found the Lekner method ^22] of image cells 
particularly straight forward. The equations of motion 
are solved by the adaptive-step fifth-order Runge-Kutta 
method |23| . The dislocation number is constant, since 
so far we considered neither dislocation creation nor an- 
nihilation. Equal numbers of dislocations with positive, 
h = -(-x, and negative, h = —x, Burgers vectors, render 
the system neutral. The dislocation collective speed (also 
called activity) V{t), is defined as, V{t) = J2iLi 
where Vi — dxi/dt. The acoustic emission signal is pro- 
portional to the dislocation collective speed. Another 
popular choice is V{t) = ^f^ibiVi{t), which is propor- 
tional to the strain rate [TT] . 

Adiabatic Increase of External Stress:- First we con- 
sider the quasi-static or adiabatic case. After randomly 
placing the dislocations in the square cell, we let the sys- 
tem relax to the nearest (metastable) equilibrium state. 
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FIG. 2: (color online) (top) Time dependence of the collective 
speed of all A'^ = 64 dislocations in a square box of side L = 
100. (bottom) The stress plotted against the total dislocation 
displacement for the same run. Displacement at time t is the 
total distance all the dislocations traveled from the beginning 
of the simulation (t — 0) till time t: dt' 'Yla^i bidxi{t'). The 
arrows indicate the start and finish of the last large avalanche. 
(Note that in the stress vs displacement bottom figure the 
equilibration occurs at zero external stress.) 



During that procedure we apply zero external stress. As 
the system approaches the nearest energy minimum the 
dislocation motion slows down. A simple eigenmode 
analysis shows that the time needed for the system to 
reach zero activity diverges. We assume that the system 
is sufficiently close to the energy minimum when the dis- 
location activity has fallen below a threshold, Vth = 0.1, 
which is roughly 100 times less than the initial activity 
of the iV = 64 dislocation system. Once the system's 
activity has fallen below the specified threshold we start 
increasing the external stress adiabatically slowly. As 
soon as the adiabatically slowly increased stress pushes 
the system's activity above the threshold, V{t) > Vth and 
the system produces an avalanche, we keep the external 
stress constant until the avalanche stops (Fig. |2]). The 
scaling behavior is insensitive to the threshold for a value 
up to ten times larger and smaller. 

The system starts with small avalanches and as the 
stress T approaches the fiow stress Tc, it responds with 
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FIG. 3: (color online) (top) The cumulative distribution of 
critical stresses Tc for 5 different numerical densities. Each 
curve is extracted from 288 runs with A*' = 64, 32, 16 disloca- 
tions and 96 runs with A'^ = 48, 24 dislocations in a square 
box of side L = 100. The smaller the density, the nar- 
rower the distribution and smaller the mean Tc (see next 
Fig. [4|. (bottom) We obtain a good collapse using the 
curves with the three larger densities based on the expres- 
sion p{tc, p) ~ P'^ f[Tcp^ ]■ The collapse quantifies the fact 
that the distributions get steeper and have a smaller mean 
for lower density p. A = since the cumulative probability 
is restricted in [0, 1]. A' = —0.5 ± 0.02. The rescaling of the 
horizontal axis indicates that Tc ~ p"'^. 



larger and larger avalanches until at it finally flows 
steadily with an infinite avalanche. For r > Tc, the dis- 
locations keep moving indefinitely, exiting from one side 
of the simulation cell and reemerging at the other due to 
the periodic boundary conditions, without ever getting 
jammed (pinned) again. In a deformation experiment, 
this is the point when the sample yields. In summary, 
for T < Tc the system is jammed (pinned). For t > Tc 
the system is constantly flowing (yielding) (Fig. [2]). 

Jamming:- The critical stress Tc is not a universal 
quantity and every system with the same number of dis- 
locations and box size has a different Tc. We performed 
an adiabatically slow increase of the stress for different 
dislocation densities, p = N jl? . The cumulative distri- 
butions of the critical stresses is shown in top Fig. |3] 
One can observe that the distributions become narrower 



for smaller densities, as does the mean critical stress of 
the ensemble. The scaling collapse shown in bottom Fig. 
|3] gives the relationship Tc ^ ^/p. 
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FIG. 4: (color online) The mean critical stress (tc) plot- 
ted against the inverse numerical density p. The system is 
jammed for r < Tc and unjammed above it. Each point is 
extracted from 288 runs with N — 64, 32, 16, 8 dislocations 
and 96 runs with A'^ — 48, 24, 12 dislocations in a square box 
of side L = 100. The critical stress has a similar qualita- 
tive dependence on the density as in the proposed jamming 
phase diagram by Liu and Nagel [T7]. However for disloca- 
tions Tc{p) > for all non-zero densities p. 

The dislocation system exhibits jamming for r < Tc 
analogous to the work of Liu and Nagel and O'Hern 
et al. [THl HH] . Their systems are different from ours in 
that they had exclusively short range interactions (con- 
tact interactions of soft spheres) and we have long-range 
(besides the core interactions that are enforcing the "no 
climb" constraint). They observed similar distributions 
of depinning stresses and a similar concave up depen- 
dence of the flow stress on the density (Fig. |4]) . However 
in contrast to their results, we neither expect nor find a 
jamming point equivalent to their jamming point J where 
Tc = 0. This means there can be no density, however 
small, that will unjam our system at zero applied ex- 
ternal stress. The reason is that dislocations have long- 
range interactions [20]. No matter how far apart they 
are, they always feel each other. 

Theoretical Calculation of Critical Stress:- Consider 
■^R AR positive and Nj^^j^ negative edge dislocations 
parallel to the 2;-axis randomly distributed on a ring of ra- 
dius R and thickness Ai? on the {x, y)-plane. The stress 
exerted at the origin is given by, 
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adapted from [21] using 

cos(6l) cos(26') and p±(r,e') : 
We express all distances in terms of /, the mean disloca- 
tion distance, i.e. p = N/L'^ = l/l'^ in d dimensions, i.e. 



Eq.([T]) where K{e) 
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X = R/l and x = r/l. For any power law r 
tion, we get tx^ax ^ fx 



interac- 
K{9) 



N 



With p^ix,e) = ^E,=i 



X,AX S{x — Xi) 



5{e-e,). For small 



ring thickness we can approximate the integral with the 
value of the integrand at X times /S.X. The average over 
the number of the dislocations to first order gives 
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with {N+) = {N-) = (N). The effect of the number 
fluctuations on the stress per ring thickness is: 
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since are independent random variables, Poisson dis- 
tributed with the same mean and variance. Assuming 
that there are N dislocations of each kind in the entire 
area L"^ where X^ = L/l >> f , their mean number in the 
ring can be expressed as {Nx.ax) ^ x^^ ' Substi- 
tuting into Eq.a we find, {(^^ f) - ^ 



Integrating over the entire region 



with Xr: 
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0(1) the closest possible distance between 
2 dislocations gives: ^ ^ i yZJ^T^Zjni 

for 2a > d. In the thermodynamic limit, X^ — >■ oo. 



this translates to the stress scaling as 
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For 2a < d, 
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thermodynamic limit doesn't exist. For parallel straight 
edge dislocations in 2 dimensions 2a = d = 2 and 
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This agrees with our numerical result in bottom Fig. [3] 
The extraction of logarithmic corrections requires much 
larger systems than the ones that can be simulated. 

Discussion:- We were able to show, using a discrete dis- 
location dynamics model, that the mean critical stress of 
an ensemble of dislocation systems with long-range in- 
teractions. Tint ~ i/r, scales with the square root of the 
dislocation density, (tc) ~ ^/p, for straight parallel edge 
dislocations. Eq. ([t]) also agrees with the Taylor harden- 
ing relation [31] and is analogous to the effective velocity 
of a point vortex in 2 dimensional hydrodynamics |25j . 
We were able to perform the analytical calculation for 
any power law interaction, Tint ~ and for arbi- 

trary d dimensions. The theoretical result agrees with 
our simulation up to logarithmic corrections which are 
difficult to measure at system sizes amenable to simula- 
tion. Our results, both numerical and theoretical, show 



that for dislocations or particles with long-range interac- 
tions there can be no jamming point at a finite density 
(only at p — 0), provided there is no screening. 
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